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Section  1 

INTRODUCTION  AND  SUMMARY 
1 . 1 Background 

The  parabolic-equation  (PE)  technique  has  proven 
to  be  a powerful  tool  in  modeling  propagation-loss  under  a 
variety  of  circumstances.  Its  fundamental  limitation  results 
from  the  parabolic  or  small-angle  approximation  employed  in 
its  derivation.  While  the  approximation  had  been  known  for 
some  time  (Leontovich  and  Fock,  1946  and  Fock,  1946),  it 
became  truly  useful  when  integrated  with  the  split-step  FFT 
algorithm  (Tappert  and  Hardin,  1973). 

Strictly  speaking,  the  approximation  is  not  small- 
angle,  but  narrow-band.  That  is,  the  solution  will  be  accurate 
for  a small  band  of  rays  or  modes  with  phase-velocities  close 
to  a selectable  reference  (e.g.,  McDaniel,  1975  or  Fitzgerald, 
1975).  Since  in  underwater-sound  applications  small  angles 
are  nearly  always  of  interest,  this  flexibility  in  selecting 
the  accurate  band  offered  very  minimal  improvements.  As 
proposed  by  Fitzgerald  (1975)  one  could  solve  a set  of  prob- 
lems, each  centered  about  different  phase  velocities,  trans- 
ferring energy  between  appropriate  bands  as  indicated  by  a 
modal  decomposition.  While  elegant,  this  solution  requires 
solving  for  the  normal-modes  at  each  new  environment  and  is 
probably  impractical  for  most  applications. 

A similar  type  of  approach  has  been  proposed  by 
Estes  and  Fain  (1977)  where  each  Fourier  component  is  effec- 
tively propagated  with  refraction  and  reflection.  A number 
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of  approximations  and  assumptions  are  made  in  arriving  at 
what  appears  to  be  a feasible  implementation.  The  accuracy 
of  the  technique  remains  to  be  evaluated,  as  well  as  its 
practicality.  It  may,  however,  prove  to  be  a useful  refer- 
ence or  standard  for  comparison  with  approximate  techniques 
even  if  it  is  found  to  be  impractical  for  routine  applications. 

Other  attempts  at  overcoming  the  small-angle  limita- 
tion have  included  direct  integration  of  the  elliptic  wave 
equation  on  very  large,  fast  computers  (e.g.,  Brock  on  NRL's 
TI-ASC,  unpublished),  and  modification  of  the  parabolic 
operator  (Tappert,  1977).  Brock's  work  is  still  in  progress 
and  may  prove  feasible  for  limited  domains  of  high  angles 
given  such  a powerful  computer.  Tappert  has  derived  an 
operator  which  will  accurately  treat  large  angles  and  weak 
gradients,  or  small  angles  and  strong  gradients,  but  not 
large  angles  and  strong  gradients.  More  precise  limits  for 
this  new  operator  have  not  been  determined. 

1 . 2 Approach  and  Summary 

The  work  reported  here  represents  an  approximate 
solution  to  this  last  problem — large  angles  in  the  presence 
of  strong  gradients — formulated  in  a way  to  utilize  the 
efficiency  of  the  split-step  FFT  algorithm  and  hence  to  be 
applicable  over  large  ranges.  It  is  similar  to  the  approxi- 
mate solution  developed  to  remove  the  bulk  of  the  parabolic 
error  for  water-borne  paths  (Brock,  Buchal,  and  Spofford, 

1977).  Their  approach  ("CMOD")  was  to  transform  the  initial 
problem  (specifically  the  sound-speed  profiles)  into  a problem 
whose  solution  using  PE  would  be  much  closer  to  the  desired 
elliptic-equation  solution  of  the  original  problem. 
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As  is  shown  subsequently,  this  transformation 
breaks  down  (i.e.,  becomes  double-valued)  when  applied  to 
the  strong  gradients  found  in  the  refracting  sediments  of 
the  ocean  bottom.  Fortunately,  the  principal  CMOD  objective — 
convergence  zone  spacings — is  not  relevant  to  the  higher- 
angle  bottom-bounce  paths.  For  these  paths  a more  appropriate 
constraint  would  be  to  preserve  ray  cycle  distances  or  periods, 
thus  ensuring  that  each  ray  would  interact  with  the  correct 
portion  of  the  bottom  and  accumulate  the  right  number  of 
interactions  at  any  given  range. 

It  is  possible  to  achieve  this  objective  with  con- 
trollable, small,  errors  so  long  as  the  bottom  interaction 
process  can  be  treated  as  refraction  through  strong-gradient 
sediments.  This  is  typically  the  case  in  areas  where  the 
bottom  is  a good  "reflector"  such  as  abyssal  plains.  It  would 
not  apply  to  steep  slopes  of  exposed  basalt  where  the  inter- 
action process  is,  in  fact,  reflection. 

tl  c 

Hence  the  specific  problem  of  interest  is  K propa- 
gation of  high  angles  over  large  distances  where  the  energy 
is  refracted  in  the  ocean-bottom  sediments  (rather  than 
reflected  from  a hard  interface).  The  principal  difficulty 
with  PE  in  such  geometries  is  that  a steep  angle  effectively 
propagates  with  the  period  of  a shallower  angle  since  the 
horizontal  component  of  its  phase  velocity  is  less  than  it 
should  be.  Hence  the  ray  period,  or  cycle  distance,  is 
too  large.  For  water-borne  paths  this  problem  was  largely 
removed  by  CMOD.  The  relaxation  of  the  CMOD  convergence- 
zone  constraint  permits  certain  liberties  to  be  taken  with 
the  sound-speed  profile  in  the  bottom  which  could  not  be 
considered  in  the  water  column. 


For  the  bottom- interact ing  paths,  a new  algorithm, 
"BMOD,"  is  proposed  whose  objective  is  to  preserve  ray  periods 
given  the  user-specified  sound-speed  structure  in  the  bottom. 
BMOD  modifies  the  sound-speed  profile  in  the  bottom  to  compen- 
sate for  the  ray-period  error  accumulated  in  the  water  column. 
As  a result,  the  impact  of  range-varying  bathymetry  and  re- 
flectivity is  apparent  on  the  correct  rays,  and  mean  levels 
may  be  accurately  modeled.  The  position  of  a ray  within  the 
water  column  may,  however,  be  in  error  at  any  one  depth,  but 
the  error  does  not  grow  with  range.  Multipath  interference 
effects  in  the  field  will  have  an  inherent  (but  predictable) 
error  since  the  phase  velocities  of  the  high-angle  paths  will 
still  be  incorrect. 
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The  work  reported  here  has  focused  on  developing 
and  testing  the  feasibility  of  this  algorithm.  Its  full-scale 
implementation  in  a general  PE  code  (with  CMOD)  was  beyond 
the  scope  of  the  present  effort.  Such  an  implementation 
should  present  no  major  problems,  and  a number  of  associated 
issues  are  discussed  in  the  final  section  of  this  report. 

In  summary,  an  algorithm  has  been  developed  which 
is  able  to  modify  ray  periods  in  such  a way  that  high-angle 
paths  may  be  propagated  in  PE  (via  the  split-step  FFT  algorithm) 
with  sufficient  accuracy  to  model  many  important  bottom- 
related  aspects  of  propagation  loss.  The  algorithm  requires 
no  modifications  to  existing  PE  codes  since  it  basically  is 
used  as  a pre-processor  of  the  environmental  inputs. 
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Section  2 

THE  ALGORITHM  BMOD 


In  this  section  the  basic  algorithm,  "BMOD,"  is 
developed.  The  previously  developed  water-column  modification, 
"CMOD,"  is  shown  to  fail  when  applied  to  the  strong  gradients 
formed  in  the  ocean  bottom.  While  it  approximately  conserves 
ray  periods  for  water-borne  paths,  this  was  not  its  original 
intent — more  emphasis  being  placed  on  mode  turning  points. 

Since  we  are  not  generally  interested  in  sources  or  receivers 
deep  in  the  bottom,  the  turning-point  constraint  can  be  removed, 
with  ray  periods  becoming  the  principal  concern. 

2 . 1 Preserving  Ray  Periods 

The  CMOD  algorithm  approximately  preserves  the  range  variation 
of  the  true  solution  of  the  wave  equation  (WE),  whereas  PE 
errors  accumulate  with  range.  This  becomes  clear  if  we  inte- 
grate the  respective  ray  equations.  Using  Eqs.  (C24)  - 
(C26)  of  Ref.  1,  for  c = c(z)  we  have 


WE  rays  ^ — c(c^-C) 


PECMOD 


rays  £ cfc‘- C*) 


First , consider  an  ocean  of  constant  sound  speed  gradient  g 
and  of  infinite  depth  extending  downward  from  z = 0.  With 
c(z)  = cQ  + gz,  g > 0,  the  ray  passing  through  z = 0 at  an 
angle  of  0 from  the  horizontal  turns  at  the  depth  z.  = c -cn/g, 
where  c^Cq/cos©  is  the  sound  speed  at  the  turning  point. 

For  WE  the.  range  half-period  is  then 


whereas  for  PECMOD  we  find 


) 


Dm  ^ 

zjz  g L - Sc^cl . 


where  s = sec0  = c^/Cq.  Then  the  difference  in  range  half- 
periods is 


e.  -/5=S  J I 


— ^ /c1  - 

9 K 


sli'i) 

l 


+ JS- 


I 


Expanding  the  In  factor  in  terms  of  e = /(c^  - c0)/ct> 

\a<  -Jf  -j=:\  + f0(*) 


(7) 


Thus  A(0)  - 0,  and  for  small  e (small  9)  it  is  clear  that 
A(9)  >_  0.  For,  consider  the  function 


-/ f - 


f'(c) 


- -X  \ n c 

ac,  y u Vc 


c>  0 


Since  f(ct)  =0,  0 < cQ  < ct  implies  f(c  ) > 0. 


Moreover,  if  we  consider  the  ratio 


(3) 


» 


» 


A‘(e)  _ / s 

AT  +ys~- 1 

/*(«*-!) 

(9) 


then  as  9 -*■  ir/2,  rc(9)  -*■  1//2  ~ 0.707  (the  last  term  vanishes 
like  (lnx)/x  as  x = 4s  + ®).  Thus,  PECMOD  underestimates 
the  WE  ray  period  by  an  increasing  amount  as  9 increases 
from  0 to  tt/2,  but  never  by  more  than  30%.  For  RR  ray  paths 
in  the  water  column,  the  maximum  error  is  much  less.  For 
example,  Eq.  (3)  gives  rc(16°)  = 0.9967,  i.e.,  an  error  of 
only  0.3%. 
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On  the  other  hand,  the  usual  PE  algorithm  over- 
estimates the  WE  ray  period,  and  by  a much  larger  amount  than 
PECMOD  underestimates  it.  From  Eq.  (3),  the  PE  ray  period  is 
simply  a factor  of  c^/Cq  times  the  WE  period;  thus, 


The  error  in  periods  is 


seed  - Aw(<9) 


(10) 


M0)55  Af(*)  - /^(e)  = y/c;-c;  (-£-  - 1)  (id 

As  c + « (9  -*■  tt / 2 ) , A -*■  00  and  the  error  is  unbounded.  For 

t p 

ct  s Cq,  the  ratio  of  PECMOD  error  to  PE  error  is  (cf.  Eq.  (6)) 

EM  ■ -IS)  • ) 

where  s = sec9  = Ct/Cg.  Thus,  as  9 + 0,  E(9)  -*■  1/4.  For 

very  small  angles,  then,  the  PECMOD  underestimate  is  about 

1/4  of  the  PE  overestimate,  but  this  fraction  rapidly 

decreases  as  9 increases.  Returning  to  our  earlier  example, 

at  9 - 16°  A(9)  = 1 . 040A  ( 9 ) and  hence  E(9)  = 0.0033/0.040  = 
p w 

0.08. 


Next,  consider  the  typical  case  of  a water  column 
given  by  a piecewise  linear  sound  speed  profile  with  small 
gradients.  The  above  analysis  can  be  generalized  to  show 
that  RR  and  RSR  rays  confined  to  the  water  column  have  a 
shorter  period  under  PECMOD  than  under  WE.  This  fact  can 
now  be  used  to  our  advantage. 
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The  approach  of  BMOD  is  to  replace  a given  ocean 
bottom  and  sound  speed  profile  by  a pseudo  bottom  composed 
of  layers  and  a modified,  piecewise-linear  profile  such  that 
the  period  of  the  ray  turning  at  the  center  depth  of  each 
layer  is  preserved.  Extending  CMOD  directly  into  the  bottom 
might  appear  to  be  an  attractive  solution.  Since  PECMOD 
underestimates  ray  periods  in  the  water  column,  this  modified 
bottom  profile  would,  at  least  in  the  first  few  layers, 
actually  have  gradients  weaker  than  the  gradients  at  the  same 
depth  in  the  given  bottom  profile. 

Unfortunately,  this  method  fails  because  of  a 
fundamental  complication,  namely  the  essential  depth  dependence 
of  the  CMOD  algorithm.  In  fact,  for  deep  bottom  layers  and 
large  positive  gradients,  the  CMOD  depth  transformation  will 
become  double-valued!  The  transformed  depth  coordinate  z is 

z — z n M (13) 

where  n(z)  = cn/c(z)  is  the  refractive  index.  For  suffi- 
u 1/2 

ciently  large  gradients,  n 1 (z)  will  decrease  faster  than 
z increases,  and  multiple  values  of  z will  be  mapped  to  the 
same  z.  Consider  a bottom  layer  beginning  with  sound  speed 
c at  depth  z and  sound  speed  gradient  g,  i.e., 
z ^ z.  If  we  solve  the  equation 

* nVt(z)  + -f-  rifc(z)^(r)=  0 

dz  2 

_ * 

we  find  that  z will  reach  a maximum  at  z = z , 

z*=  -*). 


c(z)=c+g(z-z) , 


(14) 


where 
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and  then  z will  turn  back  on  itself.  Of  course,  only  for 
sufficiently  large  g,  namely 

a)  Hr  > <16) 

* _ 

does  z > z exist,  which  is  why  this  problem  never  occurs 
in  the  water  column.  However,  it  can  easily  happen  if  we 
try  to  exten_  CMOD  into  the  bottom.  For  example,  consider 
a bottom  layer  beginning  at  z = 4,878  m with  c = 1,524  m/ 
sec.  Then  if  g > 5/8  sec-''',  which  is  often  true  of  real 
bottoms,  CMOD  fails.  Thus,  while  CMOD  is  quite  good  for  the 
water  column  (n  =:  1),  it  cannot  handle  both  steep  angles  and 
n <<  1 . 


I 
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Fortunately,  the  principal  objective  of  CMOD — 
correct  convergence-zone  spacings — is  not  relevant  to  steep 
bottom-interacting  paths.  Our  concern  is  to  preserve  ray 
periods,  and  we  can  achieve  this  by  annexing  a well-chosen 
bottom  to  the  CMOD-transformed  water  column.  The  resulting 
problem  is  then  solved  using  PE.  The  ray  period  for  a bottom- 
interacting path  is  then  found  by  integrating  the  PECMOD 
ray  equation  (2)  through  the  water  column  and  then  the  PE 
ray  equation  (3)  through  the  bottom.  For  any  given  ray,  the 
contribution  from  the  water  column  is  predetermined  by  CMOD, 
but  BMOD  then  recursively  modifies  the  sound  speed  gradient 
in  each  layer  so  that  the  period  of  the  ray  turning  at  the 
midpoint  of  the  layer  is  preserved. 

2 . 2 Formal  Description  of  the  Algorithm 

Let  the  water  column  be  composed  of  (N-l)  layers 
such  that  c(z)  = ci  at  the  top  z = zi  of  the  ith  layer  and 
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c(z)  = ci+1  at  the  bottom  z = zi+1>  for  i=l,  N-l.  The 

reference  sound  speed  is  still  denoted  by  cQ;  in  the  following 
we  assume  Cq  is  also  the  minimum  sound  speed,  so  that  is 
invariant  under  CMOD.  Let  gi  be  the  gradient  in  the  ith  layer, 
so  that  ci+1  = ci  + §i(zi+i  - z^‘  At  the  bottom  of  the  water 
column,  z = z^  and  c = cN>  Let  the  bottom  sound  speed  have 
a given  gradient  g,  so  that  c(z)  = cN  + g(z  - zN)  for  z > zN> 
(The  extension  to  a general  piecewise  linear  bottom  is  straight- 
forward. ) Now  applying  CMOD  to  the  water  column  yields  z ^ and 
c^,  i=l,  ...,  N.  We  extend  this  modified  profile  into  the 
bottom  by  defining  layer  boundaries  ZN+1.  z>j+2  ’ •••  and  grad- 
ients g^,  gN+i>  •••  such  that  the  period  of  the  ray  turning 
at  the  midpoint  of  each  layer  is  preserved  when  PE  is  used 
to  solve  this  pseudo  problem. 

The  method  of  solution  is  recursive:  first,  g^  is 
found;  then,  gN+1  is  found  given  that  steep  rays  must  pass 
through  the  layer  with  gN,  and  so  on.  Thus,  assume  g.,  i=N, 

. . . , M-l  have  already  been  found.  Let  us  reference  rays  in 
untransformed  (not  CMOD)  space  by  their  angle  0^  at  the  sound 
channel  axis  z^,  or  equivalently,  by  the  sound  speed  c^  = 
Cq/cos0q  at  their  turning  depth.  Then  Wv  and  W^,  the  WE 
contributions  to  the  ray  half-period  from  the  water  column 
and  bottom,  respectively,  are  given  by  (cf.  Eq.  (4)). 


-f-V"*  --*•<<  * . **  *•«%**#*#$>*•*• 
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To  compute  the  contribution  C to  the  ray  half-period  from 
PECMOD,  it  is  convenient  to  integrate  in  the  untransformed 
space  with  its  piecewise-linear  profile.  Just  as  in  going 
from  Eq.  (2)  to  Eq.  (5),  we  find 


(19) 


N -1 


i-'i 


■ 

yrt+yc7c~ 

,yc7-/ct-c. 


• v/Cfc  y/c+ 

JZ  + jCx-C.r  1 . 


To  compute  the  PE  contribution  P from  the  bottom 
layers  previously  defined  by  BMOD,  we  integrate  Eq.  (3)  but 
in  transformed  space.  With  c^,  N < i <_  M,  given  by 

cl-  — (20) 


we  obtain 


H-l  ~ 

O V 

■ Vj 

II 

Ql_ 

ct  - i 

(21) 


but  now  we  must  define  ct,  the  turning  point  speed  of  the 
transformed  ray.  Its  initial  angle  0Q  at  zQ  may  be  found 
from 
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dr 


(22) 


& 

dr  dr  dr  dz. 

where  dz/dr  is  obtained  from  Eq . (2),  and  dz/dz  from  Eq.  (14) 
Thus , 


A 0, 


i-f-K  — /2c  -£) 


(23) 


Then  the  turning  point  speed  ct  is  given  by  Snell's  Law  in 
the  transformed  medium: 


L<5 

ct  = — - 

6 


(24) 


Note  that  c^.  is  not  the  same  as  the  image  of  c^  under  CMOD. 

The  final  ingredient  required  to  equate  ray  periods 

t h 

is  the  contribution  from  the  current  (M  ) layer.  Assuming 
that  ct  > cM,  the  Mth  layer  with  gradient  gM  will  contribute 


T = 


c. 


3- 


A/  , 

c;-c. 


(25) 


-c  nm 


to  the  ray  half-period.  Then  g^  can  be  obtained  by  equating 
half-periods : 


V + V “ c + p + 


(26) 


or  more  explicitly, 
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(?) 


& ya1-  c.v 

c.(VJ„+Vi-C-P) 


The  problem  now  is  to  find  the  unique  0Q  and  gM 
such  that  Eq.  (27)  is  satisfied  and  the  turning  depth 


Zk 


KV 

z + 

4 — Ite 


A/ 

C»-G, 

v 

9» 


(28) 


is  also  the  midpoint  l/2(z^  + z^+^)  of  the  next  bottom  layer. 


That  is , we  want 


az(5) 


%) 

<V/ 

1 - c 

Ay 

(?) 

1 

This  transcendental  equation  for  9q  can  be  solved  iteratively 

by  the  secant  method.  Let  9^.^,  i=l,  2,  . ..,  denote  the 

u ( 1 ) 

sequence  of  iterates;  corresponding  anlges  9q  may  be  found 
by  solving  Eq.  (23)  for  c^  and  using  Snell's  Law.  We  choose 

Co 


*Jl»  ~i  / Co  \ 

*•  - **  ( —) 


(30) 


to  be  the  initial  angle  of  the  ray  which  turns  at  the  top 
of  the  layer;  then  Az(0q1^)  = 0.  Next,  choose  6q2^  to  be 
the  initial  angle  of  the  ray  which  would  turn  at  the  midpoint 
of  the  current  layer  if  the  previous  gradient  were  extended; 
i.e.  , 


e. . 


- COS*  £ 


] 


(31) 
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Then  compute 

az<  r ) - 

using  Eq.  (24)  to  obtain  ct,  and  Eq.  (27)  for  g . Now  the 
secant  method  consists  of  applying  repeatedly  the  scheme 


I 


I 


J)  - /±z[daU '* ) 


/ ^ Ci-I  ) \ A I \ 

Az(0  )-Az(ffo  ) 


a.  t ; -i) 


(33) 


for  i=3,  4,  and  with  D given  by  Eq.  (29).  The  iteration 

terminates  when 


A z (O-  D 

for  some  tolerance  factor  £ chosen  in  advance;  e = 0.01  worked 
well  in  practice,  and  typically  at  most  five  iterations  (i  <_  5) 
were  required  for  convergence. 

Finallly,  gy  is  taken  to  be  the  gradient  found  on 

the  last  iteration.  This  algorithm  could  fail  only  if 

W + W.  - C - P < 0 for  all  9n  > 9~^,  in  which  case  it 
w b 0 0 

would  be  impossible  to  find  a positive  gM  from  Eq.  (27). 

But  in  Section  2.1  we  saw  that  W + W,  - C > 0 for  all  9-.. 

_ w b 0 

Also,  as  9 q -*■  ir/2  (ct  -*■»),  P -*■  0 in  Eq.  (20).  Thus,  for 
sufficiently  large  9Q,  Ww  + - C - P > 0 and  gM(9)  > 0, 

so  there  is  no  problem. 


< 


(34) 
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As  a practical  matter,  it  is  good  that  at  the  bottom 
grazing  angle,  PECMOD  underestimates  the  ray  period.  Thus, 
there  is  no  initial  angular  gap  (beyond  the  grazing  angle) 
where  ray  periods  cannot  by  preserved.  Then  by  including 
enough  bottom  layers,  the  spacing  of  the  angles  where  periods 
are  preserved  is  narrow  enough  so  that  period  errors  at 
intermediate  angles  are  small. 

The  generalization  of  the  algorithm  to  the  case  of 
multiple  bottom  layers  (gradients)  in  the  given  ocean  is 
immediate . 
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Section  3 

PE  IMPLEMENTATION  OF  BMOD 

In  this  section  we  present  a PE  implementation  of 
BMOD  for  a homogeneous  ocean  and  show  that  good  results  were 
obtained.  We  close  with  a discussion  of  some  issues  which 
would  have  to  be  addressed  in  a full  implementation  of  BMOD 
in  production  PE  codes. 

3 . 1 Homogeneous  Ocean 

Although  a full  PE  implementation  of  BMOD  was 

beyond  the  scope  of  this  contract,  it  was  feasible  to  examine 

the  simplest  case  of  an  isovelocity  water  column  c(z)  = Cq, 

since  then  the  CMOD  transformation  reduces  to  the  identity 

map.  The  following  such  case  was  tested  numerically:  cn  = 

-1  u 

1524  m/sec,  bottom  sound  speed  gradient  g = 0.5  sec  , 

source  depth  z = 12.2  m,  water  depth  3930  m,  acoustic 
s 

frequency  25  Hz,  and  a source  aperture  of  30°. 

The  new  PE  source  function  of  Garon,  Hanna,  and 
Rost  (1977)  was  used  to  illuminate  the  desired  angular 
aperture,  beyond  which  the  intensity  fell  off  rapidly  enough 
that,  at  least  for  the  first  ray  period,  the  path  of  the 
steepest  ray  could  be  fairly  well  identified.  Thus  the  ray 
period  of  the  steepest  ray  could  be  well  estimated  from 
contour  plots  of  the  acoustic  intensity  produced  by  the  PE 
runs. 

7 

The  PE  depth  mesh  interval  was  45.7  m,  and  2 
mesh  points  were  required.  The  BMOD-modif ied  bottom  sound- 
speed  profile  is  shown  in  Table  1 at  half-mesh  intervals. 
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Table  1 

MODIFIED  BOTTOM  SOUND  SPEEDS  ( BMOD ) 


Depth 

(m) 

Sound  Speed 
( ft/sec) 

3930.3 

1524.4 

3953.2 

1556.4 

3976.0 

1588.5 

3998.9 

1637.  2 

4021.7 

1686. 0 

4044.6 

1738.  8 

4067.4 

1791.7 

«v 


I 


I 


» 


» 


rf(e.)-  & rjt)  - rjt, ) (36) 

Co 

Thus  we  find 

7^(3 0°)  « 17.30  km,  rfC3oc)  - 19.99  km 

f 

In  the  actual  PE  run  with  BMOD,  the  30°  period  appeared  to 
be  about  17.6  km,  or  much  closer  to  WE  than  to  PE. 

• 3.2  The  General  Inhomogeneous  Case 

The  numerical  testing  of  BMOD  which  fell  within  the 
scope  of  this  contract  has  indicated  its  applicability  to  a 

* more  general  class  of  inhomogeneous  problems.  As  bottom 
gradients  and  water  depths  increase,  however,  we  must  expect 
that  BMOD  will  become  less  efficient.  For  example,  as  the 
ratios  of  sound  speeds  at  adjacent  PE  mesh  points  increase 

i 


I 


The  midpoint  of  each  layer  is  identified  by  the  angle  (in 

the  water  column)  of  the  ray  which  turns  there  and  whose  period 

is  preserved.  Below  the  last  depth  shown  a layer  of  constant 

sound  speed  was  added  to  insure  no  steeper  bottom  returns. 

Now  the  WE  period  for  a ray  with  angle  0 in  the  water  column 

s 

is 

-*[  A 

where  H is  the  depth  of  the  water  column.  The  PE  period  for 
the  same  ray  (unmodified  profile)  is  then  obtained  merely  by 
a multiplication: 


+■  0* 


3 


(35) 
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(i.e.,  as  bottom  gradients  increase),  the  split-step  Fourier 
algorithm  suffers  more  and  more  from  aliasing ; that  is,  the 
fraction  of  total  energy  which  is  aliased  will  grow.  The 
only  remedy  for  this  problem  is  to  increase  the  sampling  in 
depth,  i.e.,  to  increase  the  number  of  depth  mesh  points. 
However,  this  doubly  reduces  computational  efficiency:  not 
only  do  the  FFTs  in  depth  take  longer,  but  the  range  step  Ar 
must  be  decreased. 

Another  possible  limitation  is  that  the  PE  depth 
mesh  spacing  controls  the  spacing  of  the  angles  of  the  rays 
whose  periods  are  preserved.  If  that  angular  spacing  gets 
too  wide,  significant  errors  in  period  could  occur  for 
intermediate  angles. 

If  the  modified  gradients  become  too  steep,  the 
rays  of  interest  will  want  to  turn  back  up  in  the  bottom 
almost  instantaneously.  Clearly,  a ray  which  spends  a hori- 
zontal distance  less  than  Ar  in  the  bottom  cannot  be  treated 
accurately  by  PE.  Taking  Ar  very  small  is  computationally 
inefficient . 

Finally,  it  appears  feasible  to  ixtend  BMOD  to 
slowly  varying  range-dependent  environments,  but  again,  care- 
ful numerical  testing  would  be  required. 
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Section  4 
CONCLUSION 

We  have  presented  a new  algorithm  BMOD  for  use  in 
solving  for  acoustic  pressure  fields  with  significant  energy 
returned  from  the  ocean  bottom.  When  combined  with  the  pre- 
viosly  introduced  CMOD,  BMOD  creates  a pseudo-problem  which, 
when  solved  by  PE,  should  be  a good  approximation  to  the  true 
WE  solution  for  a certain  class  of  problems.  This  class 
should  be  determined  by  a full  PE  implementation  of  BMOD, 
followed  by  numerical  evaluation. 
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